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Abstract 

This paper describes the results from a collaborative study 
between the European Space Operations Center, the European 
Organization for the Exploitation of Meteorological Satellites, the 
National Oceanic and Atmospheric Administration, and the Coop- 
erative Institute for Meteorological Satellite Studies investigating the 
relationship between satellite-derived monthly mean fields of wind 
and humidity in the upper troposphere for March 1994. Three 
geostationary meteorological satellites GOES-7, Meteosat-3, and 
Meteosat-5 are used to cover an area from roughly 1 60° W to 50° E. 
The wind fields are derived from tracking features in successive 
images of upper-tropospheric water vapor (WV) as depicted in the 
6.5-(i absorption band. The upper-tropospheric relative humidity 
(UTH) is inferred from measured water vapor radiances with a 
physical retrieval scheme based on radiative forward calculations. 

Quantitative information on large-scale circulation patterns in the 
upper troposphere is possible with the dense spatial coverage of the 
WV wind vectors. The monthly mean wind field is used to estimate 
the large-scale divergence; values range between about -5 x 10~ 6 
and 5 x 10 -6 seer 1 when averaged over a scale length of about 
1000-2000 km. The spatial patterns of the UTH field and the 
divergence of the wind field closely resemble one another, suggest- 
ing that UTH patterns are principally determined by the large-scale 
circulation. 

Since the upper-tropospheric humidity absorbs upwelling radia- 
tion from lower-tropospheric levels and therefore contributes signifi- 
cantly to the atmospheric greenhouse effect, this work implies that 
studies on the climate relevance of water vapor should include three- 
dimensional modeling of the atmospheric dynamics. The fields of 
UTH and WV winds are useful parameters for a climate-monitoring 
system based on satellite data. The results from this 1 -month 
analysis suggest the desirability of further GOES and Meteosat 
studies to characterize the changes in the upper-tropospheric mois- 
ture sources and sinks over the past decade. 
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1. Introduction 

The water vapor band of the geostationary satel- 
lites is routinely used for extracting upper-tropospheric 
humidity (UTH) and water vapor (WV) windiields. 
Apart from the utility of both products in numerical 
weather prediction and model validation, it has be- 
come evident that UTH and WV winds constitute an 
important source of climatological information. Only a 
geostationary platform provides measurements of 
water vapor amount and motion on a large scale from 
a single observing system; thus, these measurements 
sen/e as truly independent datasets for comparisons 
with other climatologies. Through joint efforts at the 
European Space Operations Center (ESOC) and the 
Cooperative Institute for Meteorological Satellite Stud- 
ies (CIMSS), data from three geostationary satellites 
(GOES-7 a t 112°W, Meteosat-3 at 75°W, and 
Meteosat-5 at the Greenwich Meridian) are processed 
for March 1994 with similar algorithms to derive UTH. 
and WV winds covering 1 60°W-50°E (Fig. 1 ). Monthly' 
mean observations are computed from two UTH and 
WV wind fields per day (1200 and 0000 UTC) for 
GOES-7 and Meteosat-3 ; for Meteosat-5 eight UTH 
and four WV wind fields are used per day. 


2. Importance of UTH 

Satellite determinations of UTH have been shown 
to be useful for deriving climatologies (van de Berg et 
al. 1991; Wu et al. 1993) that can be compared with 
those derived from models (Del Genio et al. 1994; 
Stephens et al. 1993; Soden and Bretherton 1993; 
Schmetz and van de Berg 1994). The comparisons 
reveal significant differences between UTH from sat- 
ellite observations and those in general circulation 
models. 

Figure 2 illustrates that all altitude levels of the 
water vapor profile, including those represented by the 


1578 


Vol. 76, No. 9, September 1995 



Fig. 1. Coverage from the three geostationary satellites 
(GOES-7 at 1 12°W, Meteosat-3 at 75°W, and Meteosat-5 at the 
Greenwich Meridian) out to 50° from nadir. 


UTH, are important to the water vapor greenhouse 
effect. The greenhouse effect (e.g., Inamdar and 
Ramanathan 1 994) is defined here as the difference 
. between the longwave flux emitted by the surface and 
the outgoing radiative flux (OLR) escaping from the 
top of the atmosphere. A radiative transfer model is 
used to compute the change of the OLR as a function 
of perturbations in the vertical water vapor profile. A 
10% increase in relative humidity is added to a stan- 
dard tropical profile at a single level at a time; the levels 
are at 1-, 3-, 5-, 7-, 9-, 1 1 and 13-km altitudes. The 
sensitivity is expressed as the decrease in OLR di- 
vided by the increase in total precipitable water, 
del(OLR)/del(PW). Figure 2 shows that this sensitivity 
increases from about 2.8 to about 1200 W nr 2 (g cnrr 2 )- 1 
as the 10% relative humidity perturbation is shifted 
from 1- to 13-km altitude (solid curve). A mixing ratio 
profile typical for a tropical atmosphere is also shown 
in Fig. 2 (dashed curve). The greenhouse effective- 
ness of a single H 2 0 molecule (solid curve) rapidly 
increases with altitude. This is to be expected since 
radiation from the lower atmosphere absorbed by H 2 0 
molecules high in the troposphere is reemitted at 
coldertemperatures, in contrastto reemission at warm 
temperatures by the lower-atmospheric H 2 0 molecules. 
The exponential decrease of the water vapor with 
altitude balances the increase of the greenhouse 
effectiveness of a single water vapor molecule. The 
combined effect suggests that the sensitivity of OLR to 
perturbations in the water vapor profile in a tropical 
atmosphere is roughly constant up to an altitude of 
— about -10-kmr-Thus— Figr-2-indieates-the-need-for- 
observations of the UTH, as well as observations of 
. total precipitable water (which is dominated by low- 
level water vapor) in studies of the atmospheric green- 
house effect. 


3. UTH algorithms 

a. Meteosat-5 

The UTH is inferred from the WV band (5. 7-6.9 /i), 
which has a spatial resolution of 5 km x 5 km. The WV 
band observes the earth’s atmosphere in a spectral 
region characterized by strong H 2 0 absorption and is 
sensitive to approximately the upper 3 mm of total 
water vapor (Schmetz and van de Berg 1994). In a 
cloud-free atmosphere, this WV band senses typically 
down to a pressure level of 500-600 hPa. The height 
allocation of the broad contribution function depends 
on the actual profiles of temperature and humidity and 
the satellite viewing angle (Fischer et al. 1981). 

The UTH retrieval follows the approach described 
by Schmetz and Turpeinen (1988). Radiative forward 
calculations are performed for different relative hu- 
midity profiles. The necessary temperature profiles 
are taken from short-term model forecasts provided 
by the European Centre for Medium-Range Weather 
Forecasts. The different relative humidity profiles are 
constructed as follows: below 600 hPa the forecast 
profile is used, between 600 and 300 hPa the relative 
humidity is fixed at one of six values (1%, 5%, 10%, 
20%, 40%, 100%), and above 300 hPa the relative 
humidity decreases linearly to 0% at 100 hPa. Com- 
parisons of the UTH retrieval (van de Berg et al. 1 991 ) 
with radiosondes show a root-mean-square differ- 
ence of less than 1 0%. The Meteosat-5 UTH calcula- 
tions use clear-sky radiances obtained from a multi- 
spectral image classification for segment areas 
(32 x 32 pixels), ranging from 1 60 up to 230 km on a 
side. 


dOLR / dPW [ W m J /(g cm*’)] Mixing Ratio (g/kg) 



Pig. 2. Sensitivity of the outgoing longwave radiative, flux density^ 

(OLR in W nr 2 ) to changes in precipitable water (PW in g cm* 2 ). PW 
has been varied at different levels, and the sensitivity is expressed 
as the decrease in OLR over the increase in PW. The dashed curve 
corresponds to the mixing ratio profile characteristic of a tropical 
standard atmosphere. 
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The Meteosat-5 UTH in this study is somewhat 
different from that presented in van de Berg et al. 
(1991). That work confined valid UTH retrievals to 
segment areas that were totally free of medium- and 
high-level cloud so that the UTH was representative of 
the whole segment. Here the UTH retrieval is extended 
to segments that are partially covered with high- and 
medium-level cloud by using the WV radiances from 
the portion of the segment that contains clear sky. This 
tends to moisten the UTH, especially in the tropical 
convective regions where high-level clouds prevail. 

The operational calibration of the Meteosat WV 
band is based on radiative transfer calculations with 
radiosonde profiles and collocated clear-sky WV ob- 
servations. A recent recalibration algorithm has been 
added for this study, which improves calibration per- 
formance (Schmetz and van de Berg 1994). 

b. Meteosat-3 

Meteosat-3 data are similar to those of Meteosat-5, 
differing only in spectral range (5.8-7. 3 p). However, 
the UTH retrieval from Meteosat-3 differs from the 
Meteosat-5 UTH retrieval in several aspects. The 
major difference is that the Meteosat-3 UTH is re- 
trieved for each pixel regardless of cloud conditions, 
so it has much higher spatial resolution. When the WV 
radiance for a given pixel is affected by cloud, the UTH 
does not represent a column-mean relative humidity 
and can indicate up to 1 00% relative humidity. Includ- 
ing cloudy pixels in the retrieval increases the monthly 
mean UJH by about 10% (absolute) in convective 
areas; no change is observed in areas without high- 
level cloud. As a result, monthly mean Meteosat-3 
UTH values are slightly more moist than those from 
Meteosat-5 in cloudy regions. 

There are also several minor algorithm changes. 
Soden and Bretherton (1993) showed that, under 
normal conditions, the WV band brightness tempera- 
ture is close to a logarithmic function of the UTH. 
Accordingly, the forward calculation of UTH from 
Meteosat-3, which is used to generate the lookuptable 
for WV vapor band brightness temperature, is now 
performed on a different set of values (1 %, 5%, 1 5%, 
30%, 60%, 100%) and the interpolation is changed 
from linear to loglinear. Finally, the CIMSS forward 
model is used for radiative transfer computation. 

c. GOES-7 

The GOES-7 WV band (6.4-7. l/i) has a spatial 
resolution of 7 km x 7 km at nadir. While the UTH 
retrieval algorithm used with GOES-7 is identical to 
that used with Meteosat-3, the UTH fields from the 
two satellites are somewhat different because the two 
instruments sense in a slightly different part of the 
spectrum. In a clear atmosphere, the GOES-7 WV 
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band typically senses down to a pressure level of 
400-500 hPa, depending on view angle and atmo- 
spheric moisture content (Hayden 1988). Thus, com- 
pared to the Meteosats, theGOES-7WV band senses 
radiation from a slightly higher layer in the upper 
troposphere. For March 1 994 over an area (30°N- 
30°S, 82°W-105°W) with relatively equal and small 
viewing angles for both GOES-7 and Meteosat-3, 
tests showed the GOES-7 UTH was typically 5% 
higher than the Meteosat-3 UTH. 

The empirical distribution function (EDF), described 
by Weinreb et al. (1989), is employed to eliminate this 
intersatellite bias. For each of the concurrent observa- 
tions from the above-mentioned area (62 for March), 
an EDF is constructed that is used to adjust the 
GOES-7 UTH to have the same statistical distribution 
as the MeteosatS UTH. The unadjusted GOES-7 
UTH, which is representative of a higher layer in the 
upper troposphere, is consistently more moist than the 
Meteosats UTH. This suggests that the relative hu- 
midity for the period and area under consideration 
increases with height in the upper troposphere. Al- 
though such profiles are not uncommon, it is generally 
agreed that relative humidity decreases with height for 
longer-term averages. A possible explanation is that 
Meteosats calibration was found to be warm with 
respect to GOES-7, which would make the UTH drier 
(Menzel et al. 1994). Further investigation is required 
to clarify this issue. 

d. Mean UTH for March 1994 

Figure 3 shows the monthly mean UTH for March 
1994 from all three satellites. The adjusted GOES-7 
UTH merges very well with the Meteosat-3 UTH with- 
out noticeable discontinuity, while the transition from 
Meteosats to Meteosat-S is distinguishable (prima- 
rily because different UTH retrieval procedures are 
used). A unified retrieval for all satellites would likely 
alleviate these boundaries. 

Highest relative humidities of about 60%-70% are 
observed over the tropical convective zones of South 
Africa and South America. The relative humidity val- 
ues are computed with respect to water vapor satura- 
tion pressure. (At an altitude of 400 hPa, the ratio of 
saturation pressure over ice to that over water is about 
0.8.) 


4. WV wind fields 

a. Meteosats 

The method for deriving Meteosats WV winds at 
ESOC is described by Laurent (1993) and Holmlund 
(1993). The height assignment of the WV winds fol- 
lows two approaches, (a) If the target area for tracking 
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Fig. 3. The monthly mean UTH for March 1994 from all three satellites. Values have been averaged to roughly a 1° lat-long grid. The 
adjusted GOES-7 UTH merges very well with the Meteosat-3 UTH without noticeable discontinuity. On the other hand, the merge line 
between Meteosat-3 and Meteosat-5 is distinguishable, primarily because different UTH retrieval procedures are used. 


the displacement contains cloudy features, the height 
assignment is the same as for the high-level infrared 
cloud motion winds (Schmetz et al. 1993). Opaque 
clouds are assigned to the pressure level where the 
brightness temperature fits a short-term forecast pro- 
file; semitransparent clouds are assigned to a pres- 
sure level using the IR and WV bands, (b) If the target 
area is cloud free, the coldest 25% of the WV pixels in 
a target area is used to compute a brightness tempera- 
ture for the height assignment (Laurent 1993). The 
brightness temperatures are converted to pressure 
heights via short-term forecast temperature profiles. 
The adequacy of this approach has been verified in 
detailed comparisons with rawinsondes (Holmlund 
1993). 

b. Meteosat-3 and GOES-7 

Water vapor wind sets derived at CIMSS from 
Meteosat-3 and GOES-7 are processed from hourly 
full-disk images using fully automated procedures 
[see Merrill et al. (1991) for further details]. Potential 
targets are identified by maximum brightness tem- 
perature gradients exceeding predetermined thresh- 
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olds within a selected target region. Heights are as- 
signed using the WV band brightness temperature of 
the coldest 25% of the target region matched to a 
short-term forecast profile. The vertical distribution 
of the pressure heights is typically in the range of 
200-500 hPa with a peak near 300 hPa. As a final 
step, an objective quality control procedure is invoked 
once the wind vector field has been derived (Hayden 
and Velden 1991). Statistical comparisons between 
WV wind vectors derived at CIMSS and collocated 
Western Hemisphere rawinsondes show rms differ- 
ences that are close to those of current operational 
upper-air cloud-drift winds (Velden et al. 1 994). Water 
vapor winds have also been shown to have positive 
impact on NWP (Velden 1995). 

c. Mean WV wind field for March 1994 

The daily WV winds produced by both CIMSS 
(2 d ') and ESOC (4 d ') for March 1994 are used to 
calculate a monthly mean field by averaging all wind 
vectors contained within 2.0° lat-long geographical 
regions (grid points). A monthly mean gridpoint motion 
vector value is considered to be valid when the num- 
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Fig. 4. The mean WV wind vectors displayed in kts for March 1 994. (Vectors every 4° are shown.) Blue indicates vectors with pressures 
greater than 400 hPa; brown indicates vectors with pressures between 300 and 399 hPa; and green corresponds to vectors with pressures 
less than 300 hPa. 


ber of observations exceeds 20. (Only grid points near 
the limb fail this criterion.) The monthly mean gridpoint 
height value is calculated by averaging individual 
vector height assignments. Figure 4 illustrates that the 
mean WV wind vectors show good consistency and 
the heights are confined to adequately describe the 
upper-tropospheric layer mean flow. 


5. UTH and WV wind field divergence 

To associate the observed upper-tropospheric cir- 
culation with the observed UTH, divergence fields are 
computed with a finite-difference scheme from 
the monthly mean WV winds. The dashed contours in 
Fig. 5 indicate mean divergent areas that are often 
associated with water vapor emerging from source 
regions (deep convective areas). For example, the 
pattern in the UTH over South America is clearly 
correlated with the divergence field. 


The solid contours in Fig. 5 correspond to upper- 
tropospheric convergence, most prominent over sub- 
tropical subsidence regions. The observed convergence 
values appear quite reasonable: assuming that the 
convergence takes place in a layer of about 1 50 hPa, 
the associated subsidence is about 60 hPa d~\ which 
compares favorably to the observed subsidence in the 
trade wind regions (Betts and Ridgway 1989). 

A comment is warranted concerning the height 
variation of the WV winds. Inspection of Fig. 4 shows 
that the monthly mean altitudes of the WV winds range 
approximately between 250 and 350 hPa. The spatial 
changes are smooth, and the standard deviation at a 
geographical location is only 30-70 hPa. Thus, the 
computed divergence/convergence patterns do not 
correspond to a constant pressure surface but to a 
smoothly varying surface with higher altitudes in the 
convective area than in the subsidence regions. 

One should not expect a perfect match between 
UTH and wind field divergence. Over tropical areas 
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Fig. 5. Simultaneous display of the upper-level convergence fields (derived from WV winds) and the observed UTH for March 1 994. The 
dashed contours correspond to divergent areas in the upper troposphere; solid contours represent convergence. Contours are shown in 
intervals of 2 x 1 0' 6 s 1 . The UTH is displayed as a color image with a relative humidity scale displayed. (Blue is dry, green is in the middle, 
and red is moist.) 


especially, the UTH maxima are associated with con- 
vection usually forced by surface heating and/or con- 
vergence. On the other hand, the WV wind field is 
primarily representative of upper-level large-scale flow, 
with secondary influence from large convective re- 
gimes. The UTH and divergence maxima should be in 
close proximity, but they need not be aligned. An 
example is evident in Fig. 5 over Africa, where the UTH 
maximum is shifted to the south of the wind divergence 
maximum. The position of the divergence maximum 
appears to be influenced by the entrance region of the 
subtropical jet over Africa. However, over the South 
American continent, where no jet stream is observed, 
the divergence maximum matches well with the UTH 
maximum. 

The congruence of the UTH and WV wind diver- 
gence patterns suggests a physical relationship in 
which the UTH patterns are largely influenced by 
upper-tropospheric transport. Picon and Desbois 
(1990) indicate that the upper-tropospheric humidity 
field is chiefly determined by large-scale dynamics. 
The clear-sky UTH climatology of van de Berg et al. 
(1 991 ) corroborates this statement, as the geographi- 
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cal distribution of monthly UTH fields is reminiscent of 
the Hadley and Walker circulation patterns. Soden 
and Bretherton (1993) also point out that the lower 
troposphere is the prime source of upper-tropospheric 
water vapor and that deep convection is the transport 
mechanism. However, the novel message of this 
study is that both UTH and the upper-tropospheric 
dynamics can be inferred simultaneously from geo- 
stationary satellites. As such, the satellite-derived 
UTH and WV winds constitute a valuable source for 
climatological diagnostics and for comparison with 
large-scale circulation models. 


6. Summary 

Monthly mean fields of the upper-tropospheric rela- 
tive humidity (UTH) and water vapor (WV) winds from 
Meteosat-3, Meteosat-5, and GOES-7 WV band ob- 
servations have been analyzed for March 1994. The 
fields provide a concise description of monthly mean 
upper-tropospheric moisture and motion fields. The 
scope of the present study is a qualitative demonstra- 
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tion that UTH and upper-tropospheric dynamics are 
directly observable on a nearly global scale; a future 
quantitative analysis will include a much larger dataset. 

The comparison of the March 1994 UTH and WV 
wind fields shows the influence of atmospheric dy- 
namics on the distribution of the upper-tropospheric 
moisture. Since upper-tropospheric moisture is an 
important factor in the atmospheric greenhouse ef- 
fect, the satellite-observed fields of UTH and WV wind 
are important for climate model validation. Most of the 
previous comparisons of climate models with ob- 
served moisture fields used only total precipitable 
water vapor observations over the oceans. Since low- 
level moisture dominates precipitable water, such 
comparisons do not validate the upper-tropospheric 
. moisture fields in models. 

This study using satellite measurements docu- 
ments an example where the upper-tropospheric 
moisture distribution is strongly influenced by the 
large-scale dynamics. Longer-term datasets of both 
UTH and WV winds covering seasonal and interannual 
variations would be useful for comparison with climate 
models and for validation of the dynamical aspects of 
the models. Such a study could help reveal how the 
large-scale circulation would be affected by climate 
change. 

This study also suggests that use of WV winds 
should improve representation of upper-tropospheric 
moisture fields in numerical models. This may be 
important to medium-range forecasts since the radia- 
tive heating/cooling critically depends on upper-tropo- 
spheric moisture. The importance is amplified by the 
fact that the moisture field influences the cloud field, 
which has an even larger effect on radiative heating/ 
cooling and also implies energy transformations 
through phase changes. 

Both UTH and WV winds could be very desirable 
contributors to a climate-monitoring system. The moni- 
toring of both quantities will soon be possible on a 
nearly global basis since future geostationary satel- 
lites from China and Japan will also be capable of 
radiance measurements in the WV band. Further 
detailed studies are needed to fully understand the 
features in simultaneous UTH and WV wind fields. A 
continuous operational derivation of the two products 
would provide the basis for investigation of seasonal 
and interannual variations. 
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